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Abstract 

We present two pieces of evidence in support of the conjecture that the 
microscopic spectral density of the Dirac operator is a universal quan¬ 
tity. First, we compare lattice data to predictions from random matrix 
theory. Second, we show that the functional form of the microscopic 
spectral correlations remains unchanged in random matrix models which 
take account of finite temperature. Furthermore, we present a random 
matrix model for the chiral phase transition in which all Matsubara 
frequencies are included. 

1 Introduction 

Motivated by the Leutwyler-Smilga sum rules for the inverse powers of the 
eigenvalues of the Dirac operator in a hnite volume 0], Shuryak and Ver- 
baarschot P introduced the so-called microscopic spectral ps density of the 
Dirac operator ip. This quantity is dehned as 

where p(A) = ~ ^n) is the spectral density of the Dirac operator, V is 

the space-time volume, and S is the chiral condensate. The microscopic spec¬ 
tral density carries information on how the thermodynamic limit is approached 
which, among other things, is very useful for the analysis of lattice data, in 
particular in the chiral limit. It was conjectured that ps should be a universal 
quantity, i.e., that it is insensitive to the details of the dynamics. As a uni¬ 
versal quantity, ps should depend only on the symmetries of the problem and, 
therefore, be calculable in a much simpler theory than QGD, namely random 




matrix theory (RMT). In RMT, the matrix of the Dirac operator (including a 
mass term) in a chiral basis in Euclidean space is replaced by a random matrix 


according to 


im Up' 


’ im W' 

.{ilpp im_ 


_W^ im _ 


( 2 ) 


where W is a random matrix of dimension N. The average over gauge held con- 
hgurations is replaced by the average over random matrices, and the weighting 
function exp(—S'giue) is replaced by a simple Gaussian distribution P(W) ~ 
exp(—Tr WW"^) of the random matrices. Depending on the representation 
of the fermions and the number of colors, there are three different universality 
classes which have been classihed in Ref. [^. They are described by the three 
chiral Gaussian ensembles: the chiral Gaussian Orthogonal Ensemble (chGOE) 
where W is real, the chiral Gaussian Unitary Ensemble (chGUE) where W is 
complex, and the chiral Gaussian Symplectic Ensemble (chGSE) where W is 
quaternion real. The microscopic spectral density has been computed for all 
three ensembles in the framework of RMT iS- 

There were a number of Endings which added support to the universality 
conjecture for ps. We refer the reader to a recent review by Verbaarschot 
[^] which also gives a comprehensive account of other applications of RMT 
in QGD. In this paper, we present direct evidence for the universality of ps 
by comparing lattice data to RMT predictions. This is done in Sec. pT| . In 
Sec. we add another piece of evidence by showing that the functional form 
of ps and of the microscopic spectral correlations in general does not change 
in random matrix models which take account of finite temperature. In Sec. 
we discuss a model for the finite-temperature chiral phase transition in which 
all Matsubara frequencies are included. Gonclusions are drawn in Sec. |^. 


2 Universality of and of the microscopic spectral cor¬ 
relations 

2.1 Lattice data and predictions from random matrix theory 

We analyze lattice data obtained by Berbenni-Bitsch and Meyer for an SU(2) 
gauge theory with staggered fermions in the quenched approximation with 
(5 = 2.0. To study finite-volume effects, data were obtained for four different 
lattice sizes: 4^, 6^, and 10^. For each lattice, the complete spectrum 
of the Dirac operator was obtained for a very large number of configurations 
(9979, 9953, 3896, and 477 configurations for the 4^, 6^, 8^, and 10^ lattice, 
respectively). This was necessary to obtain good statistics. The numerical 






eigenvalues of the Dirac operator satisfy an analytical sum rule with excellent 
precision (relative accuracy ~ 10“®). 

According to Ref. [Q, staggered fermions in SU(2) should be described 
by the chGSE. Analytical results exist for the distribution of the smallest 
eigenvalue (only for the quenched approximation) |]^, 

^(Amin) = y|c(cA„,i„)®/^/ 3 / 2 (cAmi„)e“^(‘=^-")' , (3) 

and for the microscopic spectral density (also for the unquenched case) [|], 

1 1 

ps{z) = 27r(c^)^y duu^ J dv [Ja-l{2uvcz)ja{2ucz) — vja-l{2ucz)ja{2uvcz)] . 
0 0 

(4) 

Hence, a direct comparison with lattice data is possible. Here, a = iVj + 2z/ + l, 
where iV/ is the number of massless dynamical quarks and u is the topological 
charge. The parameter c sets the scale on which the eigenvalue is measured. 
It should be pointed out that c is not a free parameter but determined by the 
lattice volume and the chiral condensate according to Eq. ©• 

In Figure 1 we have plotted the data and the predictions from RMT for 
both P(Amin) and Ps{z) for all available lattice sizes except for the 10^ lattice 
where the number of conhgurations was too small to obtain sufficiently good 
statistics. The parameter a in (^) is 1 since Nf = 0 and since = 0 on 
the lattice [§]. The data agree almost perfectly with the RMT predictions for 
all but the smallest lattice. The range over which ps agrees with the RMT 
prediction increases with the lattice size which is evident by the dehnition of 
Ps in Eq. (|^). We conclude that a lattice size of 6"^ is sufficient to identify the 
universal microscopic spectral density of the Dirac operator. 

Unquenched calculations and their analysis in terms of RMT are in pro¬ 
gress. It will be very interesting to see the effect of a small dynamical quark 
mass in pg. 

2.2 The microscopic spectral density at finite temperature 

Various random matrix models have been constructed for the spectrum of the 
Dirac operator at hnite temperature to describe generic features of the chiral 
phase transition. We shall discuss these models in more detail in Sec. The 
common feature of all these models is that a temperature-dependent matrix Y 
is added to the random matrix W in Eq. (^. The specihc form of Y depends 
on the choice of basis. However, as we shall see in Sec. ^ the matrix Y can 
always be chosen real and diagonal. This dehnes the generic problem. 
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Figure 1: Distribution of the smallest eigenvalue P{Xram) and microscopic 
spectral density Ps{z) for three different lattice sizes (see text). Note that 
there are no free parameters in the RMT prediction. 








The matrix of the Dirac operator (for m = 0) thus has the form 


0 W + V 
+ Y 0 


( 5 ) 


with Y = diag(Yi,..., Yjy). The are real and depend on T. The dimension 
N of the matrix W can be identihed with the space-time volume V. For the 
moment, let us concentrate on the chGUE and on the quenched approximation. 
In this case, we can compute the chiral condensate which is given by 

EiT) = J:^x, ( 6 ) 


where S is the chiral condensate at T = 0 and x is the only real and positive 
solution of 

1 ^ 1 

= (7) 


N + y? 


or zero if no such solution exists 0. The computation of the microscopic spec¬ 
tral density and the corresponding microscopic spectral correlations is more 
difficult since the standard methods involving orthogonal polynomials cannot 
be applied here due to the presence of the additional matrix Y. We have used 


the graded eigenvalue method [|T3 and extended it to the present case with the 
additional chiral symmetry of the Dirac operator MM- A similar calculation 
has been done simultaneously and independently by the authors of Ref. 


la 


The result of the calculation is very simple: The functional form of ps and 
of the microscopic spectral correlations remains unchanged provided that we 
take into account the temperature dependence of the chiral condensate 
Specihcally, we hnd that the fc-point functions can be written as 


14 


Rk{xi, ...,Xk)= det[G(a;p, Xg)]p^g=i^,„^k , 


( 8 ) 


where in the microscopic limit the function C is given by 


C{xp,Xg) = 2N 


^ UpTi (Up)T q(U g) UqjQi^Up') Ji(^Uq') 


lU 


ut 


u„ 


(9) 


with Up = 2NExp and Ug = 2NExq. Hence, the only change with respect 
to the zero-temperature case is that the microscopic variables are rescaled by 
S, the temperature-dependent generalization of S. This result adds further 
support to the conjecture that ps (and the microscopic spectral correlations) 
are universal. 





















So far, the calculation has been done for the quenched approximation only. 
The unquenched case gives rise to some technical complications, but we have no 
doubt that the same result will be obtained: The functional form of the micro¬ 
scopic spectral density and the microscopic spectral correlations will not change 
provided that the arguments are rescaled by the temperature-dependent chiral 
condensate. We make the same conjecture for the chGOE and the chGSE 
although it is not clear at the moment how the calculation should be done in 
these two cases since the graded eigenvalue method requires the computation of 
an Itzykson-Zuber-like integral which is not yet known for the chGOE and the 
chGSE. Since the microscopic spectral correlations are known at zero temper¬ 
ature for all three chiral ensembles (even in the presence of massless dynamical 
quarks) and since the temperature dependence of the chiral condensate can be 
computed easily for each ensemble, the hnite-temperature results thus follow 
immediately. 

3 A random matrix model for the chiral phase transi¬ 
tion with all Matsubara frequencies included 

As we have already mentioned in the previous section, at hnite temperature 
a temperature-dependent matrix Y has to be added to the random matrix W 
representing the matrix elements of the Dirac operator. Why this is so will 
become obvious below. Since P{W) is invariant under a transformation of 
the basis, we can always make Y diagonal by a suitable basis transformation. 
While at zero temperature the choice of basis was completely irrelevant we 
now have to specify the basis with respect to which the matrix elements of the 
Dirac operator are evaluated. A different choice of basis will be reflected in a 
different form of the matrix Y. In Refs. |T5| , a plane wave basis was chosen 
whereas the basis of Ref. @ consisted of (anti-) instanton zero modes. In 
only the lowest Matsubara frequency was included in Y since it is sufficient 
to determine the properties of the chiral phase transition at T = T^. We now 
present a model which includes all Matsubara frequencies, thereby extending 
the validity of the model to lower temperatures. [| (For a comparison with the 
model presented in Ref. |]^ see below.) 

If one chooses a plane-wave basis 0fc(x)e*'^''’' at hnite temperature, the do- 
term in ip gives rise to a term proportional to ujk = {‘2k + l)7rT {k integer) in 
the matrix of the Dirac operator. These terms are included in Y. The basis 

^After completion of this work, we were informed by M. A. Nowak that an equivalent 
model, derived from the NJL model and additionally including a chemical potential, was 
discussed recently in Ref. 0 - One can show that the two models yield identical results. 
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states 0fc(x) corresponding to the same Matsubara frequency uik are coupled 
by a random matrix. In the same way, basis states belonging to different 
Matsubara frequency should be coupled by random matrices. The matrix of 
the Dirac operator at hnite temperature can thus be written in the form 

im hh + Dl , . 

.lyt + D im J ^ 

with 


D = diag( 7 rT,..., vrT, —ttT, ..., —ttT, SttT, ..., SvrT, —3ttT, ..., —SttT, ...) . 

( 11 ) 

The dimension of a submatrix of D belonging to one single Matsubara fre¬ 
quency is equal to K which should be identihed with the three-volume V 3 . The 
four-volume is V 4 = V 3 /T = K/T. The elements of the random matrix W are 
distributed according to P{W) ~ exp(—Tr hTITl). The model presented 
in Ref. also took into account all Matsubara frequencies. However, basis 


states corresponding to different Matsubara frequencies were coupled by the 
same random matrix. There is no reason why this should be the case, and the 
resulting temperature dependence of the chiral condensate shows non-analytic 
features which are physically unrealistic. 

Going through the same steps as in |^, we obtain a saddle-point equation 
similar to (|^ from which the chiral condensate follows as (qq) = This 

equation reads (in the chiral limit m = 0 ) 




1 


N 


E 

n=l 



2K 1 

K/T to 


( 12 ) 


where the upper limit N/2K of the sum should go to inhnity. We now perform 
the sum, choose units in which {qq){T = 0) = 1, and express the result in terms 
of {qq) to obtain 


{qq) = tanh 


(99) 
T/T, • 


(13) 


This equation can only be solved numerically, and the result is plotted in Fig¬ 
ure 2. The low-temperature behavior of the chiral condensate can be obtained 
analytically in this model. We have 


(qq) ^1 — 2 exp 



as T —> 0 . 


(14) 


We observe that the chiral condensate does not change signihcantly unless the 
temperature is close to T^. While this is in qualitative agreement with lattice 












Figure 2: The chiral condensate in units of its value at T = 0 for the model 
discussed in the text. 


data 1^1 Eq. (H) does not agree with the well-established low-temperature 
expansion of chiral perturbation theory []^. The reason for this discrep¬ 
ancy is that the random-matrix model in its present form is effectively zero¬ 
dimensional in space and, therefore, unable to describe the effect of the Gold- 
stone bosons adequately. In order to account for the Goldstone bosons, addi¬ 
tional non-random parts would have to be added to the random matrix. 

At the critical temperature, we can compute the critical exponents /3, 7 , 
and 5. These are defined using the reduced temperature t = {T — Tc)/Tc by 
{qq) ~ \t\^, X ~ 1^1”^ with x = d{qq)/dm\m=o, and {qq) ~ at T = Tc. We 
obtain /? = 1/2, 7 = 1, and 5 = 3 in agreement with This supports the 


position taken in that only the lowest Matsubara frequency is responsible 
for the critical behavior at the phase transition. 


4 Conclusions 

We have presented very direct evidence for the universality of the microscopic 
spectral density of the Dirac operator by comparing lattice data obtained by 
Berbenni-Bitsch and Meyer to predictions from random matrix theory. Im¬ 
pressive agreement was found already for relatively small lattices. Glearly, in 
the microscopic limit the lattice data are described by random matrix theory. 
Thus, RMT is a very useful tool in the analysis of lattice data and may provide 




















guidance for further lattice calculations. It remains to be seen whether RMT 
can also prove useful with respect to simplifying the generation of gauge held 
conhgurations. This would be enormously helpful, in particular close to the 
chiral limit. 

The universality of the microscopic spectral correlations of the Dirac op¬ 
erator was further supported by a generic hnite-temperature calculation in 
which it was found that the functional form of the microscopic spectral cor¬ 
relations does not change with temperature provided that the temperature 
dependence of the chiral condensate is taken into account in the microscopic 
scale. Although this calculation was only done for the chGUE in the quenched 
approximation, it is more than likely that analogous results will be obtained 
in the unquenched case and for the other two chiral ensembles. 

A model was presented for the chiral phase transition at finite temperature. 
This model took account of all Matsubara frequencies. It reproduced the ob¬ 
servation from lattice data that the chiral condensate remains almost constant 
up to temperatures close to and also conhrmed the position of Ref. 


that the lowest Matsubara frequency is responsible for the properties of the 
chiral phase transition. Reasons for discrepancies with the low-temperature 
expansion of chiral perturbation theory have been discussed. 

A number of very interesting topics such as the effect of a chemical potential 
[p^ , pT| and the question of topological charge and the axial U(l) anomaly 
have not even been mentioned in this article. Although the field of applying 
RMT to problems of QCD is a very recent one, many useful applications have 
already been found, and we are convinced that this will not cease to be the 
case in future work. 
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